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We study the electronic band structure and optical conductivity of an array of tunnel-coupled 
array of graphene nanoribbons. We show that due to the coupling of electronic edge states for the 
zigzag nanoribbon structure, the Fermi surface can become a hyperbola similarly to the case of 
the layered metal-dielectric structures, where the hyperbolic isofrequency contours originate from 
the coupling of localized surface plasmon polaritons. Moreover, we show that for both types of the 
ribbon edge, the optical response of the structure can be characterized by a uniaxial conductivity 
tensor, having principal components of the different signs. Therefore, the tunnel-coupled nanoribbon 
array can be regarded as a tunable hyperbolic metasurface. 


I. INTRODUCTION 

Graphene is a two-dimensional carbon lattice that ex¬ 
hibits a wide range of unique electronic properties such 
as linear dispersion of charge carriers, Berry phase, and 
room-temperature quantum Hall effect^. Beyond the 
purely electronic properties, graphene has attracted a 
lot of attention in the field of nanophotonics, since it 
is both transparent and exhibits strong light-matter in¬ 
teraction^. Moreover, graphene supports the surface 
electromagnetic waves, analogous to surface plasmon po¬ 
laritons in metallic and semiconductor films. The ad¬ 
vantage of the graphene plasmon-polaritons over their 
metallic and semiconductor counterparts is that they can 
be efficiently controlled with the external gate voltage^. 
The attractive plasmonic properties of graphene nanos¬ 
tructures led to the emerging of the field of graphene 
nanoplasmonics^. 

One of the topics in graphene plasmonics is the study 
of graphene metasurfaces, the arrays of graphene nanois¬ 
lands which may exhibit the optical properties which 
differ significantly from those of the two-dimensional 
graphene sheets^—. For example, it has been shown 
that a two-dimensional array of graphene nanoislands can 
play the role of a perfect absorber for electromagnetic 
wavesi^. In the studies of graphene-based metasurfaces 
it is commonly assumed that the conductivity of indi¬ 
vidual graphene elements coincides with the conductiv¬ 
ity of two-dimensional graphene. However, it is known 
from numerous studies that the graphene nanopattern- 
ing can substantially alter the electronic band structure 
of graphen o 16 ’ 17 . Specifically, in Ref it was shown that 
one can open and efficiently control the width of the band 
gap by nanopatterning graphene. Altering the electronic 
band structure should lead to the modification of the AC 
conductivity and thus to the modification of the optical 
properties of the nanopatterned graphene. In Ref^ it 
was theoretically shown that the quantum confinement 
effects start play the substantial role in the optical re¬ 
sponse of the patterned graphene if the characteristic size 
of the nanostructure elements becomes smaller than 10 
nm. Technically, the introduction of the superlattice to 
the graphene sheet at these scales can be realized by ei¬ 



FIG. 1. (Color online) (a) Geometry of the structure (b,c) 
Structure the array of tunnel coupled array of zigzag (b) or 
armchair (b) graphene nanoribbons. 


ther etching the graphene sheets, introducing the peri¬ 
odic strain^ - — or by introducing the periodic gate volt¬ 
ag e 23 ’ 24 . Moreover, in Ref^ the possibility of molecular 
assembly of the arrays of several-atoms-thick nanorib¬ 
bons was shown. 

In our work we are studying the structures shown in 
Fig.©, which are essentially the arrays of tunnel-coupled 
graphene nanoribbons. We demonstrate that the optical 
properties of these structures can be tailored by engi¬ 
neering the electronic band structure, and that such a 
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system can be regarded as a uniaxial anisotropic medium 
both for electrons and plasmon-polaritons. The optical 
properties of the graphene nanoribbon arrays have been 
studied in a number of papers. Specifically, in2£ it has 
been shown that a nanoribbon array can play a role of 
an efficient polarizers for the THz radiation due to the 
huge absorption anisotropy. However, in all of the pa¬ 
pers the tunneling between the individual ribbons was 
not accounted for. The allowing for the electron hop¬ 
ping between the ribbons results in a number of physical 
effects such as emergence of hyperbolic electronic bands 
from the edge states in the zigzag ribbon case and over¬ 
lap of electron and hole bands in the armchair edge case 
which will be discussed further in detail. 

The remainder of the paper is organized as follows. 
In section EH we show the results on the electronic band 
structure and Fermi surface topology of the arrays of 
zigzag and armchair graphene nanoribbons. Section Mil 
presents the results on the conductivity spectra obtained 
with the Kubo formula. The conclusion of the paper is 
given in section ITVl Finally, appendices CY1 and iBl provide 
the detailed formalism for the calculation of the band 
structure and the conductivity, respectively. 


II. BAND STRUCTURE CALCULATION 

For both types of the ribbons edge the Hamiltonian 
can be written as2I: 

h =—t y {tfjbp + b\a a ) -1 y 

nearest next 

nearest 

-T Y (aafys+^Oa), (!) 

inter— 

ribbon 

where the first term corresponds to the hopping of the 
electron between the nearest atoms with the character¬ 
istic hopping energy t = 2.8 eV, the second term to 
the next-nearest hopping with t' = 0.2 eV, and the last 
term corresponds to the tunneling between the nearest 
ribbons with the tunneling energy r which can be var¬ 
ied (see FiglTJ). Generally, r should decrease exponen¬ 
tially with distance between the neighbouring ribbons a 
t ss t exp(—2(d — a)/a)2£, where a is the interatomic dis¬ 
tance. We note that as a goes to infinity, the band struc¬ 
ture and the conductivity should reduce to those of an 
individual nanoribbon. In the opposite limit a = a, the 
considered structure can either converge to the conven¬ 
tional 2d graphene case as in Fig[Tfb) or result in a more 
complicated structure different from the two-dimensional 
graphene as in FigOIc). In our calculations we have used 
a = y/3a which corresponds to t = t'. 

Indices [ a , /?] in the Eq. (JT]) each consist of three inte¬ 
gers: the first one i is the number of the atom in the unit 
cell. The unit cells for the case of armchair and zigzag 
ribbons are shown with dashed lines in Figs.[]Jb,c). Each 
unit cell contains N atoms of each sublattice. The re¬ 
maining two indices (m, n ) correspond to the position of 
the unit cell in the directions of K b and k , respectively 


(cf. Figs [Ub, c) for further details). We then introduce 
Fourier transform of the annihilation and creation oper¬ 
ators a a ,a) a : 

= 4, a i {k,K b )e~ i ^ mL+kna \ 

N c k,K b 

«Un = E 4 K b )e i( ' KbmL+kna ' > (2) 

N <= k,K b 


where N c is the number of unit cells. The expressions for 
the operators b are identical. The wave functions of the 
structure are the linear combinations of the single atom 
eigenfunctions: 


<t>(k, K b ) = Y [^i a ! (k, K b )\0) + Bib\{k, iG,)|0}] , (3) 


i=l,JV 


and thus can be represented as a 2 N vector of the expan¬ 
sion coefficients A,B: 

<P=(A 1 ...A n ,B 1 ...B n ) t . (4) 

If we substitute the expressions © to the Hamilto¬ 
nian © we obtain H = Y^ H(k,K b ), where H(k,K b ) 

k,K b 

can be represented as a 2 N x 2 N matrix in the basis of 
the single atom eigenfunctions: 


Explicit form of the matrix depends on the type of the 
ribbon edge and is presented in appendix[X]for both types 
of the edge, here we just note that due to the requirement 
of Hermiticity the condition H B a = h\ b holds. 

The problem is then reduced to the eigenvalue prob¬ 
lem for a 2 N x 2 N matrix H(k , K b )(j) = e(k, K b )</>, which 
has been solved numerically in order to obtain the band 
structures and Fermi surfaces. For the case of zigzag 
ribbons N = 10 and for the case of armchair ribbons 
N = 11. The results are presented in Figs l2l3l We first 
consider the zigzag edge case shown in Fig|2] We first 
note,that due to the inclusion of the second nearest neigh¬ 
bours hopping, the electron-hole symmetry is broken and 
the electron and hole bands touch at the nonzero en¬ 
ergy^. It can be seen that the individual branches of the 
single ribbon band diagram evolve to the bands. More¬ 
over, the degeneracy of the two edge states is lifted and 
we observe two distinct bands formed by the tunnel cou¬ 
pled edge states. These bands resemble the plasmonic 
bands existing in the layered metal-dielectric structures^ 
and originating from the photon tunneling between indi¬ 
vidual surface plasmon polaritons propagating along the 
metal-dielectric interfaces. It is these plasmonic bands 
that form the hyperbolic isofrequency contours and make 
the layered metal dielectric structures the most conven¬ 
tional realization of hyperbolic medium^. The coupled 
waveguide arrays have been show to exhibit hyperbolic 
media properties both for the case of electromagnetic ) 32 ! 33 
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FIG. 2. (Color online) Electronic band structure of the ar¬ 
ray of tunnel coupled zigzag nanoribbons, (b) Contour plots 
of the Fermi surfaces of the structure. The numbers at the 
curves define the corresponding electron energy. 


and acoustio^ waves. As can be seen in Fig(2](b), where 
the energy contour plots are presented, the similar situ¬ 
ation takes place for the electrons in coupled zigzag rib¬ 
bons. At the energy where the electron and hole bands 
touch (equal to 3 1' rts 0.6eV) the contour plots reduce to 
single points at the touching position. Then, as we shift 
to the area where the band formed by the edge states 
exists, we can observe the hyperbolic-like isofrequency 
contours. Thus, if the chemical potential is tuned such 
that the Fermi energy falls to the narrow region of these 
bands, the electrons at the Fermi surface would propa¬ 
gate in the effective hyperbolic media. As long as ef¬ 
fectively ballistic propagation of electrons is considered, 
the negative refraction and the partial focusing of elec¬ 
tron waves should be observed at the interface of the 2d 
graphene and patterned graphene structure just as for the 
case of electromagnetic waves— The hyperbolic isofre¬ 
quency contours also lead to the enhanced density of elec¬ 
tronic states as compared to the case of two-dimensional 
graphene at the same energy. We also note that at certain 
energies the Fermi contour reduces to parallel lines, and 
the two-dimensional structure effectively behaves as a one 
dimensional. As it is known^, in the vicinity of such 
points the electron-electron interactions are substantially 
modified and the inclusion of the Flubbard repulsion to 
our tight-binding model is the subject of the future work. 


The advantage of the graphene-based structures is that 
the Fermi level can be fine-tuned with the external gate 
voltage, such that it falls to the region characterized by 
the hyperbolic isofrequency contour. 

The band diagram and the Fermi surface contour plots 
for the case of an armchair ribbons are shown in Fig|3] 
For the case of armchair ribbon, the Dirac cones exist 
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FIG. 3. (Color online) Electronic band structure of the array 
of tunnel coupled armchair nanoribbons, (b) Contour plots of 
the Fermi surfaces of the structure. Blue (red) curves corre¬ 
spond to electron (hole) branches. The numbers on the curves 
define the corresponding electron energy. 

only if the number of atoms in the unit cell N = 3m + 2, 
where m is an integer—. In the considered structure N = 
11 and thus the Dirac cones are preserved. It can be 
seen that coupling leads to the overlap of the electron 
and hole bands, and thus as can be seen in Figj3](b) at 
certain Fermi energies, the Fermi surface contains both 
electron and hole pockets. 

We now move to the calculation of the AC conductivity 
of the considered structures. 


III. CONDUCTIVITY TENSOR. 


When the eigenvalues and eigenvectors are obtained, it 
is straightforward to compute the AC conductivity tensor 
of the structure using the Kubo formalism. The conduc¬ 
tivity tensor is given by22 


H (/[ £ m] “ f{ £ n}) (0ml - e«a|0 

k,K b 


(Tm - 4)K - e s n + tuo + iS) 


( 6 ) 
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The velocity operator v is defined through — ev/c = 
dH(A)/dA, where H(A) is derived from the tight- 
binding model via Peierls substitution: 

. , . In' Adr _ , 0 j^(R'-R)A ( 7 \ 

tR.R' ^R.R’ ehc R ~tRM' eKc ■ {<) 

The latter equality is made within the approximation 
that the vector potential is constant at the scale of single 
unit cell. The velocity operator then can be presented as 
an 2 N x 2N matrix of the form: 


—ev = —e 


(vaa vab 
\v B a v B b 


( 8 ) 


Due to the requirement of the hermiticity of the veloc¬ 
ity operator, we set vab = and within the nearest 
neighbour approximation vaa = vbb = 0. The explicit 
form of vab depends on the type of the edge and is 
presented in the appendix E The results are shown in 
Figs0|a,b). The temperature was set to 300K and the 
chemical potential was equal to 3 1'. The conductivity is 
normalized to the bare graphene conductivity <To = e 2 //i 
For the case of zigzag ribbons, we observe the strong ab- 



FIG. 4. (Color online) Real and imaginary parts of the prin¬ 
cipal components of the conductivity tensor for the case of 
zigzag (a) and armchair (b) ribbon edge. Letters in (a) 
and (b) correspond to the transitions shown in Fig[2](a) and 
FigEKa) respectively. 

sorption for the field polarized along the ribbons. This 
effect is due to the presence of the edge states and is also 
present in the conductivity spectra of individual zigzag 
ribbons^ and graphene nanodiscs^. However, due to the 
lift of degeneracy between the edge states, an additional 
absorption channel between the two edge states bands ap¬ 
pears, which results in the broad resonance of the absorp¬ 
tion for the transverse-polarized electric field, depicted 
with letter A in FigQJa). In the low frequency region, 


the conductivity has the Drude-like shape for both elec¬ 
tric field polarizations. In the higher frequency region 
the conductivity is determined by the interband transi¬ 
tions shown in Fig[2](a). Namely, the transitions between 
the second hole band and first electron band and first 
electron and second hole band result in the two distinct 
absorption peaks in xx conductivity tensor component 
(marked with letters C,D) and one peak in zz compo¬ 
nent (marked with B). Transitions between the first hole 
and third electron and first electron and third hole bands 
result in the higher frequency absorption peak for both 
conductivity components (marked with letters E,F). 

For the case of armchair ribbons shown in Fig. |U(b) 
it is worth noting that the interband absorption is al¬ 
most fully suppressed for the transitions between the 
first electron and hole bands. This is quite different to 
the case of conventional graphene which is characterized 
by a broadband absorption caused by interband transi¬ 
tions. The suppression of interband absorption caused 
by the Fermi surface topology has been studied previ¬ 
ously^. The detailed study of the interband absorption 
suppression in the considered geometry is the subject of 
the ongoing work. In the higher frequency region the ab¬ 
sorption peaks are defined by the interband absorption 
between the higher electron and hole bands. We note 
that for the xx component of the conductivity the peak 
in absorption is a step-like similar to the two-dimensional 
systems, and zz conductivity has the peaks similar to the 
Van Hove singularities in one-dimensional systems. For 
the electromagnetic applications, it is important that in 
the certain frequency region the absorption is suppressed, 
and the imaginary parts of the conductivity tensor com¬ 
ponents have opposite signs, which suggests that we can 
realize a high quality anisotropic metasurface based on 
the considered structure. 

We should also note, that the obtained conductivity 
tensors are nonlocal, i.e. the conductivity depends on 
the in-plane momenta of the electromagnetic wave . In 
Fig. [5] the real part of the conductivity of both zigzag and 
armchair ribbons is presented for different values of Kf,D 
and compared to the conductivity of two-dimensional 
graphene. The substantial modification of the real part 
of the conductivity is observed due to the opening of 
new absorption channels via the indirect interband tran¬ 
sitions. 


IV. CONCLUSION AND OUTLOOK. 


In this paper, we have analyzed the electronic band 
structure and conductivity spectrum of two-dimensional 
arrays of tunnel-coupled arrays of graphene nanoribbons. 
We have shown that both depend drastically on the type 
of the ribbon edges. Namely, for the case of zigzag 
edges the coupling of the electronic edge states lifts the 
degeneracy between the symmetric and antisymmetric 
electronic eigenmode formed by the edge states. More¬ 
over, each of the individual modes evolves to a narrow 
band, defined by the hyperbolic isofrequency surfaces. 
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FIG. 5. (Color online) Conductivity spectrum for the zigzag 
(a,b) and armchair (c) ribbons for different values of K b D. 
The conductivity of two-dimensional graphene is shown for 
comparison. 


We have thus demonstrated the simple yet interesting 
analogy between the considered structures and multi¬ 
layered metal-dielectric metamaterials, where the hyper¬ 
bolic isofrequency surface for the electromagnetic field 
originates from the coupling between individual surface- 
plasmon polariton modes. For the case of armchair ri- 
bons, we have shown that the coupling leads to the over¬ 
lap of electron and hole bands, resulting into the forma¬ 
tion of the hole and electron pockets in the Fermi surface 
of these structures in the certain doping range. 

Furthermore, we have shown that these structures 
can be regarded as metasurfaces for the electromagnetic 
waves in the optical to mid-IR frequency ranged. Specif¬ 
ically, we have demonstrated that both types of the 
structures are characterized by different signs of imag¬ 
inary part of conductivity tensor components together 
with small levels in some frequecy ranges. Recently, a 
phenomenological study of the electromagnetic spectrum 
of such metasurfaces has been performed, showing that 
these structures can support a special type of Dyakonov- 
like hybrid surface electromagnetic waves^. 

The modification of the electronic band structures as 
well as electronic eigenfunction profiles in the nanostruc- 
tured graphene should lead to the substantial modifica¬ 
tion of both electron-electron and electron-phonon inter¬ 


actions. Tailoring the electron-phonon scattering via the 
electron band engineering may lead to the substantial 
modification of the plasmon scattering rates and thus to 
the modification of the optical response of these struc¬ 
tures^ Self-consistent calculation of the Coulomb self¬ 
energy as well as interactions-induced modification of the 
AC conductivity is the subject of future work. Avail¬ 
ability to control the optical response of the graphene 
metasurfaces via the modification of the electronic band 
structure through the quantum confinement makes these 
structure a perspective platform for the novel optical to 
IR optoelectronic devices. 


Appendix A: Band structure calculation 

Below is presented the explicit form of the Hamiltonian 
from Eq. © for both types of the ribbons edge. 

Zigzag edge: 


Hab = —t ■ 


/k 0 0 
1 K* 0 
0 1k 


0 ^ e iKbL \ 
0 0 

0 0 


(Al) 


Haa 


0 0 

0 ... 

K 0 


Vo 0 

0 ... 

1 K* 
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/ cos ka 
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2 
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.. 0 \ 
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cos ka 

Kj* 

2 

0 

0 
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2 

cos ka . 

0 

0 

0 

0 
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2 

V o 

0 

0 

.. cos ka J 


(A2) 


where k = 1 + e ika , L = NaV^/2 H B a = H\ b ,H bb = 

zr* 

h aa- 

Armchair edge: 


Haa = 
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-If- 



/I 
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0 0\ 


1 e 

ky/3a 


0 0 

—t ■ 
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Hbb = 


/ 
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r_ e iK b L\ 
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K 
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0 
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0 

0 
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where k = 1 + e ikaV ^,L = Na/ 2, H B a = h\ b . 


(A5) 


Appendix B: Calculation of the AC conductivity 


. The expression for the conductivity operator for the 
case of zigzag edge is given by: 

J 


(r 2 + r 3 e 
ri 
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ika 


ev ab = — ~rt ■ 

n 


0 

r 3 + r 2 e 
r± 


— ika 


0 
0 

r 2 + r 3 e 


ika 


given by 






0 
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0 
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0 
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a/V 3( 

-1,0); 

f2 = 

a/V 3(1/2, >/3/2) 

r 3 


1 

with 
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